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ABSTRACT 

A  numerical  model  for  heat  transfer  in  laminar  duct  flows  has  been  developed 
using  the  finite  difference  method  to  explore  the  significance  and  extent  of  "back- 
conduction"  at  low  Peclet  numbers.  The  calculations  have  been  carried  out  for  flows 
between  parallel  plates  and  in  circular  tubes  by  using  different  Peclet  numbers  in  the 
range  of  0.05  to  100.  For  both  situations  constant  heat  flux  and  constant  wall  temperature 
boundary  conditions  were  used.  The  validity  of  the  results  has  been  checked  by 
comparison  with  some  existing  results  in  the  literature,  and  extended  to  a  wider  range  of 
parameters  including  conjugate  wall  conduction  effects.  The  results  are  presented  for  bulk 
mean  temperature  variation,  Nusselt  number  behavior,  and  energy  absorbed  before  the 
heated  section,  for  cases  with  and  without  wall  conduction.  Such  axial  conduction  effects 
may  be  an  important  feature  in  the  thermal  characterization  of  microtubes,  which  are  to 
be  used  in  microheat  exchangers. 
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I.      INTRODUCTION 


The  heat  transfer  problem  of  laminar  fluid  flow  in  ducts  which  is  known  as  the 
Graetz  problem  has  many  applications  in  technology  and  has  been  studied  extensively 
since  Graetz  (1885).  The  classical  Graetz  problem  considers  the  forced  convection  heat 
transfer  of  the  fluids  flowing  in  ducts  neglecting  axial  conduction  effects  in  the  fluid  or 
the  wall. 

This  study  is  on  the  axial  conduction  effects  of  the  flows  whose  Peclet  numbers 
range  is  from  those  of  liquid  metals  where  axial  conduction  may  not  be  neglected,  to 
those  of  oils  where  the  axial  conduction  has  almost  no  effect  on  the  temperature 
distribution  inside  the  ducts.  Developing  micromachines  nowadays  and  in  the  future,  the 
axial  conduction  effects  may  be  a  very  important  feature  in  microheat  exchangers. 
The  effects  of  axial  conduction  on  the  Hagen-Poiseuille  flows  between  parallel  plates  and 
in  circular  cylinders  for  constant  wall  temperature  and  constant  heat  flux  have  been 
studied  in  this  thesis.  The  governing  equations  have  been  solved  using  the  finite 
difference  numerical  scheme.  The  validity  of  the  results  has  been  checked  by  comparison 
with  some  existing  results  in  the  literature  and  extended  to  a  wider  range  of  parameters 
including  conjugate  wall  conduction  effects.  The  results  are  presented  for  bulk  mean 
temperature  variation,  Nusselt  number  behavior,  and  energy  absorbed  before  the  heated 
section,  for  cases  with  and  without  wall  conduction. 


II.   BACKGROUND 
A.  HEAT  TRANSFER  IN  LAMINAR  DUCT  FLOWS 

Consider  a  flow  in  a  circular  tube  where  fluid  enters  the  tube  with  a  uniform 
velocity.  The  viscous  effects  are  important  in  this  flow  and  a  boundary  layer  develops  as 
x  increases.  The  flow  is  "fully  developed"  in  the  region  where  the  velocity  gradient 
doesn't  change  anymore  with  increasing  x  (Figure  1). 
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Figure  1.  Laminar  fully  developed  flow 


The  Reynolds  Number  for  the  duct  flows  is  defined  as 
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where  V  is  the  mean  fluid  velocity,  v  is  kinematic  viscosity,  v  =  — ,  and     DH   is  the 
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The  Prandtl  number  relates  the  temperature  distribution  to  the  velocity  distribution 
and  defined  as 


a 


(2.2) 


where  a  is  the  thermal  difrusivity  .  If  the  Prandtl  number  is  one,  the  velocity  and  the 
temperature  profiles  develop  together  and  at  the  same  rate.  The  Prandtl  number  is: 

Pr  «  1       for  the  liquid  metals 

Pr  =  1       for  the  gases 

Pr»l       for  the  oils 

For  the  liquid  metals,  where  the  Pr«l,  the  energy  diffusion  rate  is  much  more 
than  the  momentum  diffusion  rate.  It  is  opposite  for  the  oils,  in  which  Pr»l  and  the 
velocity  profile  develops  faster  than  the  temperature  profile  in  this  case.  The  value  of  the 
Prandtl  number  strongly  affects  the  relative  growth  of  the  velocity  and  thermal  boundary 
layers. 

Peclet  number  is  defined  as 

Pe=Re.Pr  (2.3) 

and  it  is  a  measure  of  the  quantity  of  the  axial  heat  conduction  effects  in  the  fluid.  The 
axial  conduction  is  assumed  negligible  for  Pe>10  and  the  axial  conduction  term  can  be 
assumed  small  in  the  governing  equation  for  this  case.  But  the  axial  conduction  effects 
can  be  significant  when  the  Peclet  number  is  smaller.  The  purpose  of  this  numerical  study 
is  to  show  the  importance  of  the  axial  conduction  effects  and  the  heat  absorption  upto  the 
entrance  from  where  the  heating  starts  when  the  Peclet  number  is  small. 

The  Nusselt  number  is  defined  as 

Nu  = *.  (2.4) 

k 

where  k  is  the  thermal  conductivity  of  the  fluid.  It  is  a  non-dimensional  number  where  the 
convection  coefficient  h  is  calculated. 

The  laminar  duct  flow  being  considered  is  assumed  to  be  hydrodynamically  fully 
developed  before  any  heating  effects  are  considered.  The  velocity  profile  in  the  laminar 
fully  developed  region  is: 
For  the  circular  tube: 


u(r)  =  2V 


1- 
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(2.5) 


For  the  parallel  plates: 


u(y)  =  -V 


1- 


(2.6) 


where  r  and  y  are  measured  from  the  centerline,  r0  is  tube  radius  and  h  is  the  distance 
from  the  centerline  to  the  plates.  [Ref.l] 

1.  The  Heat  Transfer 

The  heat  transfer  for  a  duct  flow  can  be  expressed  using  Newton's  law  of  cooling, 
q  =h(Tw-Tm),  where  h  is  the  convection  heat  transfer  coefficient,  Tw  is  the  wall 
temperature  and  Tm  is  the  mean  temperature  of  the  fluid,  where  Tm  is 
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For  the  constant  cv  and  incompressible  flow  through  the  circular  tube,  Tm  is 
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Vr 


JuTrdr 
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o    0 


where  for  the  flow  between  the  parallel  plates  is 
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fuTdy 
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(2.9) 


y  and  r  are  measured  from  the  centerline,  r0  is  the  radius  of  the  tube  and  h  is  the  distance 
from  the  centerline  to  the  plate.  [Ref.  1] 

The  energy  equations  for  the  cylindrical  and  rectangular  coordinates  are 
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=  0       (2.11) 


where  "i"  is  the  enthalpy.  These  equations  are  applied  to  the  problem  including  the 
boundary  conditions  and  the  temperature  distribution  is  calculated  analytically  or 
numerically.  [Ref.  2] 

2.  The  Thermally  Fully  Developed  Flow 

In  thermally  fully  developed  case,  the  relative  shape  of  the  profile  no  longer 
changes  although  the  temperature  profile  changes  with  x.  The  criterion  for  a  fully 
developed  temperature  profile  for  a  circular  tube  is 


_a_ 

dr 
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=  const  (2.12) 

r=ib 

This  condition  is  reached  in  a  duct  flow  whether  the  case  is  uniform  heat  flux  (q 
through  the  wall  is  constant)  or  uniform  wall  temperature  (The  wall  temperature  is 
constant).  These  two  cases  have  a  lot  of  applications  in  engineering.  For  example, 
constant  wall  temperature  in  boiling  or  condensation,  or  constant  heat  flux  using  the 
electrical  heater.  In  the  thermally  fully  developed  flow,  the  convection  coefficient  "h"  is 
constant,  independent  of  x.  [Ref.  1] 

B.  BACKGROUND  STUDIES 

The  heat  transfer  solution  for  laminar  fully  developed  parabolic  velocity  profile 
flow  inside  a  circular  tube  in  the  thermal  region  and  subject  to  uniform  tube  wall 
temperature  was  treated  for  the  first  time  by  Graetz  in  1883  and  is  known  in  the  heat 
transfer  literature  as  Graetz  problem.  In  his  problem  Graetz  neglected  the  axial 
conduction  and  solved  the  problem.  [Ref.  3] 

The  problem  statement  is: 

M_^)»A.£A+±«*  (2.13) 
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where 
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Tin -To 
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RenPr 


and  To  is  the  wall  temperature,  T^  is  the  initial  temperature,  r  is  the  radial  distance,  x  is 
the  axial  distance,  D  is  the  diameter,  ReD  is  the  Reynolds  number  and  Pr  is  the  Prandtl 
number.  [Ref.  4] 

The  velocity  profile  is  fully  developed  much  before  the  temperature  profile  when 
the  Prandtl  number  is  high  relative  to  one.  For  such  situations,  Graetz  solution  is 
reasonably  well  justified.  But  the  Prandtl  number  is  very  small  when  the  fluid  is  liquid 
metal  and  therefore  the  axial  conduction  effects  may  not  be  neglected.  For  such  cases, 
since  the  axial  conduction  effects  are  not  included  to  the  result,  the  solution  may  give 
inaccurate  results. 

Michelsen  and  Villadsen  investigated  the  Graetz  problem  with  axial  heat 
conduction  for  circular  tube  constant  wall  temperature  case  by  using  a  numerical 
procedure.  They  used  a  method  that  is  the  combination  of  orthogonal  collocation  and 
matrix  diagonalization.  They  didn't  include  the  wall  conduction  effects  in  the  problem. 

The  problem  may  be  defined  as: 


n       2^d0      1   a  (    80}       1     820 


dy     x  dx\    dxj     Pe    dy 


(2.14) 


where 


Pe 


2R(vz)pCi 


r 
x  =  — 

R 


y  = 


PeR 


0  =  1^ 


Tb~T0 


0  =  1^ 


Tb~T0 


The  boundary  conditions  are: 
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for     y  >  0 


where  Vz  is  the  average  fluid  velocity,  R  is  the  rube  radius,  z  is  the  axial  distance,  r  is  the 
radial  distance,  Tb  is  the  fluid  temperature  at  z  — >  -co ,  r  is  the  bulk  mean  temperature 
and  To  is  the  wall  temperature  at  z  >  0 .  They  plotted  the  results,  Nusselt  number  vs.  y, 
heat  flux  vs.  y  and  9  vs.  y  where  the  axial  heat  conduction  effects  are  easily  seen.  But 
they  didn't  include  the  axial  heat  transfer  inside  the  wall,  which  increases  the  fluid 
temperature  very  much  upto  the  entrance  as  the  wall  conductivity  gets  larger.  [Ref.  5] 

X.  Yin  and  H.  H.  Bau  included  the  axial  conduction  effects  of  the  wall  to  the 
internal  flow  through  circular  tube,  by  using  two  parameters,  duct's  outer/inner  radius 
ratio  and  fluid/wall  thermal  conductivity  ratio.  They  plotted  the  graphs  of  ""temperature 
vs.  radius  of  the  circular  tube""  and  "Nusselt  number  vs.  Peclet  number",  including  the 
axial  fluid  and  wall  conduction  effects.  [Ref.  6] 

C.  NUMERICAL  METHOD 

The  governing  equations  are  steady  state  elliptic  partial  differential  equations.  The 
temperature  T(x,y)  throughout  the  domain  must  satisfy  the  governing  equation  and  the 
boundary  conditions  along  the  entire  boundary.  The  finite  difference  method  was  used  in 
this  numerical  study  to  express  the  governing  equations  numerically.  The  central,  forward 
and  backward  difference  expressions  used  in  this  numerical  study  for  the  grid  points  in 
figure(2)  are  as  follows  : 


<- 

K 

• 

• 

• 

•j+l 

ki-U) 

• 

huj) 

• 

• 

•  J 

• 

• 

• 

• 

•j-l 

i-2 

i-1 

i 

i+1 

i+2 

Figure  (2).  The  grid  points  in  finite  difference  methods 


The  central-difference  expressions  for  the  first  and  second  derivatives  with  error 
order  of  h"  are: 


dt  _  t{i+\j)     tg-ij) 


dx              2hx 

dt  _  t(!j+l)  -  t{iJ_x) 
dy            2hy 

d  t  _  h'+ij) +  k'-ij) ' 

-2t(U) 

dx2  '               (hxf 

32 


dy2 "  {Kf 


2 


The  forward-difference  expression  with  error  order  of  h  is 

dy~  2hy 


The  backward-difference  expression  with  error  order  of  h  is: 


y 


(2.15a) 


(2.15b) 


(2.16a) 


(2.16b) 


(2.17) 


di=   hij+2)  +  *t«j+i)    *o.j)m       [Re£7]         (2>18) 

dy  2h 


Substituting  these  finite-difference  expressions  into  the  governing  equations,  the  problem 
is  reduced  to  a  set  of  linear  algebraic  equations,  depending  on  the  number  of  grid  points 
used.  Usually  large  number  of  grid  points  is  desirable,  but  the  number  of  equations  to  be 
solved  becomes  too  large  in  this  case.  Gauss  Seidel  iteration  method  can  be  used  to  obtain 


a  solution.  All  the  temperatures  of  the  grid  points  are  then  make  equal  to  the  initial 
temperature  at  the  beginning  of  the  numerical  solution.  The  temperatures  at  the  grid 
points  are  calculated  by  using  the  governing  equations  and  the  boundary  conditions.  After 
the  temperature  of  a  grid  point  is  calculated  using  the  finite  difference  equations,  this 
updated  temperature  is  used  for  the  next  grid  and  the  temperature  values  converge  to  the 
exact  temperature  a  little  more  in  every  iteration.  A  technique  called  over  relaxation  is 
used  to  speed  the  convergence  of  the  Gauss  Seidel  method  when  applied  to  the  elliptic 
partial  differential  equations.  This  technique  is: 

Hu.new)  =  Whi,j_ne*)  +  0  "  W)'(U_«W)  (219) 

w  is  known  as  the  relaxation  parameter.  For  overrelaxation  w  is  between  1  and  2;  for 
underrelaxation  between  0  and  1 .  For  this  method,  overrelaxation  must  be  used  to  speed 
up  the  convergence.  [Ref.  7] 

Forsythe  and  Wasow  show  that  for  a  45  by  45  mesh  grid  the  optimum  value  of  w 
is  around  1.870.  They  also  point  out  that  when  this  relaxation  parameter  is  used,  the 
convergence  is  approximately  30  times  faster  than  the  usual  process  (w=l).  [Ref.  8] 


10 


III.    GOVERNING  EQUATIONS 


A.  DERIVATION  OF  THE  EQUATIONS 


1.  Flow  Inside  the  Circular  Tube,  Constant  Wall  Temperature 


Case 


For  flow  in  a  circular  tube,  the  viscous  energy  equation  in  cylindrical  coordinate 


system  is 


di  di 

ox  or 


r  dr  \      dr 


\ 


1    d     .  dt 
r2  d<f>\     d(f> 


dx\    dx) 


(3.1) 


where  di  =  cdt  +  —dP  and  vr  is  the  radial  velocity  [Ref.l].  Then  the  equation  becomes 
P 

for  the  constant  k 


dt  dt 

upc v  vrpc 

dx  dr 


k  d  f   dt}     k  d2t      ,  d2t 


r  dr\    dr 


.2  cs±2 


+  k- 


rL  d<j>L        dx' 


=  0 


(3.2) 


The     temperature     distribution     is     symmetric     \d2t/d<f>2)=0     and     the     flow 


is 


hydrodynamically  fully  developed,  vr  =  0 .  Equation  then  becomes 


]_d_ 
r  dr 


f   dt} 
dr  j 


+ 


dh_ 
dx2 


=  u 


pc 


dt_ 
dx 


(3.3) 


a  =  k/pc ,  the  thermal  diffusivity  of  the  fluid.  So  the  governing  equation  for  the  circular 
duct  flow  is 


dlt     1  dt     dlt  _u  dt 
dr2     r  dr     dx2      a  dx 


(3.4) 


Now,  to  non-dimensionalize  the  problem,  assume  that 


11 


e  =  ^L 


to-te 


r+  = 


u+=-,     x+=^-,     u+=2(l-r+2) 
V  RePr 


The  governing  equation  becomes 


Aee     e2e     1  ee      1    d2e 

\-r     \——  = -  + +■ 


dx+      dr+1      r+  dr+      Pe2  dx' 


(3.5) 


2.  Flow  Inside  the  Circular  Tube,  Constant  Heat  Flux  Case 


The  governing  equation  for  the  circular  duct  flow  is  (3.4) 


d2t     1  dt     d2t      u  dt 


■  + +  ■ 


dr       r  dr     dx       a  dx 


Now,  to  non-dimensionalize  the  problem,  assume  that 


o.  '<-' 


q0d/k 


+       2xld  +      ..,       .2 

,      r '  =  — ,      u    =— ,      x    =  -         -,      u    =  2(1  -  r     ) 


V 


Re.Pr 


The  governing  equation  becomes 


+AdO       d26       1    dO        1     d26 
\-r     ]——  = -  + +  ■ 


dx+      dr+       r+  dr+      Pe2  dx' 


(3.6) 


3.  Flow  Between  the  Parallel  Plates,  Constant  Wall  Temperature 


Case 


The  viscous  energy  equation  in  cartesian  coordinate  system  for  the  flow  between 
the  parallel  plates  is 

d 


di          di           di 
up \-vp \-wp 

dx         dy  dz 


dx\    dx) 


+ 


dy 


( \  dA      d  ( ,  dt\ 

+  — 
dz 


k  — 


k 

V     dzy 


=  0      (3.7) 
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where  di  =  cdt  +  —  dP  and  vr  is  the  radial  velocity  [Ref.  1  ] .  Then  the  equation  becomes 
P 

for  the  constant  k 


dt  dt  dt      , 

upc  —  +  vpc 1-  wpc k 

dx  dy  dz 


d(dt\      d 

—   + 


dx 


'dt^ 


dyidy 


+ 


\^  J 


d_ 

dz 


(dt\ 

\dzj 


=  0 


(3.8) 


The  temperature  distribution  on  z  direction  is  constant,  —  =  0  and  the  flow  is  laminar  and 

dz 

hydrodynamically  fully  developed,  v  =  0 ,  w  =  0 .  Equation  then  becomes 


d2t     d2t 

dx2     dy2 


=  u 


( pc\dt 


V  «  J 


dx 


(3.9) 


a  -  kj pc ,  the  thermal  diffusivity  of  the  fluid.  So  the  governing  equation  for  the  parallel 
plates  flow  is 


d2t      d2t 


+ 


dxz     dy' 


u  dt 
a  dx 


(3.10) 


Now,  to  non-dimensionalize  the  problem,  assume  that 


e  =  ^- 


to-te 


+     y         +     u 
y    =—,      u    =  — ,      x 
h  V 


+  =  2^^      „+=2(iV) 

RePr  2 


The  governing  equation  becomes 


Pe 


d2e    Ad2e 

— — +  4 


u+  de 


Dh  dx+  dy+ 


2       2  dx+ 


(3.11) 
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4.  Flow  Between  the  Parallel  Plates,  Constant  Heat  Flux  Case 

The  governing  equation  for  the  parallel  plates  flow  is  (3.10) 


d2t      d2t      u  dt 


+ 


dx       dy       a  dx 
Now,  to  non-dimensionalize  the  problem,  assume  that 

q0Dh/k  h'  V'  RePr  2V     y     J 

The  governing  equation  becomes 

i     d2e     .  d2e    u+  de 

+  4 T  =  ^riTT  (3-12) 


Peh  dx+2        dy+2        2  ^ 
5.  Effects  of  Wall  Conduction,  Flow  Inside  the  Circular  Tube, 
Constant  Wall  Temperature 

The  governing  equation  for  the  circular  duct  flow  is  (3.4) 

d2t     \_dt_     d2t  _u  dt 
dr2     r  dr     dx2      a  dx 

Now,  to  non-dimensionalize  the  problem,  assume  that 

0  =  AzL,      r+=^-,      u+=^,      x+=^±,      «+=2(l-,+2) 
t„-t,  r„  V  RePr 

The  governing  equation  becomes 


AdO       d26       1    d6         1     d29 
\-r     J— —  = r  + +  — - 

'dx+      dr+~      r+  dr+      Pe1  dx+ 


(3.13) 
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B.  BOUNDARY  CONDITIONS 


1.  Flow  Inside  the  Circular  Tube,  Constant  Wall  Temperature 


Case 


Insulated  wall 

//////                 Twall=T0  =  COnStant 

x-y+<x> 
*      0->O 

_  — ._ 1 — ^. 

u(r  ) 

X 

'      '      '      '      '       '             Constant  wall  temperature 
Insulated  wall 

Figure  3.  The  boundary  conditions  of  circular  tube,  constant  wall  temperature  flow 


The  governing  equation  is  (3.5) 


M 


+2\d0     d2e     i  de      1    e2e 

\\-r     [——  = -  + + 


dx+      fr+2      r+  dr+      Pe2  dx+2 


where  the  parameters  of  the  equation  are  defined  above.  The  boundary  conditions  for  this 


case  is: 


at  r+=0, 


80 
8r  + 


=  0 


(Symmetry  boundary  condition) 


(3.14a) 


atr+=l, 


0  =  0 
80 


8r  + 


=  0 


for       x+  >  0 
f°r       x+  <  0 


(3.14b) 


as  x+->-<»,      t^>te,       0  ->  l 


as*+-»+oo,      t->tQ,       <9->0 


(3.14c) 
(3.14d) 
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2.  Flow  Inside  the  Circular  Tube,  Constant  Heat  Flux  Case 


Insulated  wall 


////// 


fl   I    1    1    fl 


'  =  0 


..........  ::±r:r: 


"u7rT 


Constant  beat  flux 
ft B fl fl fl 


e  ->  -2x" 


Figure  4.  The  boundary  conditions  of  circular  tube,  constant  heat  flux  flow 


The  governing  equation  is  (3.6) 


M 


+Ad0       d26       1    dO        1     d26 
\\-r     )— — -  = -  + +  ■ 


2         .  +    a..+  r»_2    ^    +2 


a*+       ar+        r+  5r+      /V  dx+ 


where  the  parameters  of  the  equation  are  defined  above.  The  boundary  conditions  for  this 


case  is: 


d$ 


at  r+=0,  =  0  (Symmetry  boundary  condition) 

dr  + 


(3.15a) 


atr+=l, 


86 


3r  + 
30 


ldr  + 


=  0 


1 

2 

for 

x+  >0 

for 

x+  <0 

(3.15b) 


asx+-»-oo,     t->te,      0-»O 


(3.15c) 


as  x+->+oo,      #  — »  -2x+  + 


7      r 
+ 


4  2 

+  + 


48       8         2 


(3.15d) 
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3.  Flow  Between  the  Parallel  Plates,  Constant  Wall  Temperature 


Case 


y  Insulated  wall 
/  /     /     /     /     /            Twji  =  To  =  constant 

*-*-*   „                —♦           /t    „ 

9-*\         *                   ,    u(/) 
/ 

'      '      '      '      '     '              Constant  wall  temperature 
Insulated  wall 

Figure  5.  The  boundary  conditions  of  parallel  plates,  constant  wall  temperature  flow 


The  governing  equation  is  (3.1 1) 


1      d26      A  d29      u+  d6 

+  4- 


p°lhdx+2      ay*2      2&+ 


where  the  parameters  of  the  equation  are  defined  above.  The  boundary  conditions  for  this 


case  is: 


at  y+=0, 


30 


=  0  (Symmetry  boundary  condition) 


(3.16a) 


aty+=l, 


0  =  0 
30 


3y+ 


=  0 


for       x+  >  0 
f°r       x+  <0 


(3.16b) 


as  x+  ->  -oo ,      t  -» t. ,       0  -» 1 


as*+->+oo,      /->/0,       <9->0 


(3.16c) 
(3.16d) 
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4.  Flow  Between  the  Parallel  Plates,  Constant  Heat  Flux  Case 


Z  Insulated  wall 


/ 


x->-co 
0->O 


y*  ▲      q  =<k  =  constant 


»{y) 


Constant  heat  flux 

fi 75 R ff 


Insulated  wall 


e^>-2x' 


Figure  6.  The  boundary  conditions  of  parallel  plates,  constant  heat  flux  flow 


The  governing  equation  is  (3.12) 


M 


+2\5^       d2#       1    86        1     d26> 


&+      ar+       r+  dr+      /V  &+ 


where  the  parameters  of  the  equation  are  defined  above.  The  boundary  conditions  for  this 


case  is: 


S/3 

at  y+=0,  —  =  0  (Symmetry  boundary  condition)  (3.17a) 

dr  + 


aty+=l, 


80         1 

dr+         4 

for 

x+  >0 

^_-n 

for 

x+  <0 

5r  + 


(3.17b) 


asx+-»-oo,      t->te,       <9-»0 


as  x+  ->+<»,      t^>t0,       6  =  -2x+  + 


39        3 


1120     16 


(3.17c) 
(3.17d) 
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5.  Effects  of  Wall  Conduction,  Flow  Inside  the  Circular  Tube, 


Constant  Wall  Temperature 


Insulated  wall 


Constant  wall  temperature 


//////  ^-onsiani  wan  temperature 


_cl^_ 


/'/////  7  '■    -      'P' i 

'      '      x      x      x      x  T„,,  =  TD  =  constant 


Insulated  wall 


e->o 


Figure  7.  The  boundary  conditions  of  circular  tube,  constant  wall  temperature  including  the  circular  wall 
The  governing  equation  is  (3.5) 


M 


+2\d0       d2d       1    d6        1     320 

l-r+     r— —  = -  + + 


dx+      dr+1      r+  dr+      Pe2  dx 


The  boundary  conditions  for  this  case  is: 


at  r+=0, 


30 
dr  + 


=  0 


(Symmetry  boundary  condition) 


(3.18a) 


atr+=l, 


Qyvall  ~  ® fluid 
deWall  _  d6 ) fluid 


{  dr+         dr+ 


(3.18b) 


at  r+=  !+<*•, 


d6 


dr+ 
6>  =  0 


=  0 


for      x+  <  0 
f°r      x+  >  0 


(3.18c) 


where  <5r  is  — ,  r    is  the  thickness  of  the  wall,  and  ro  is  the  inner  radius  of  the  tube. 


as  j+  ->-oo,      /->/e>       0->\ 


asx+-»+oo,      /-Mo »       0->O 


(3.18d) 
(3.18e) 
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IV.      NUMERICAL  REPRESENTATION 
A.  FLOW  INSIDE  THE  CIRCULAR  TUBE,  CONSTANT  WALL 
TEMPERATURE  CASE 

1.  The  Governing  Equation 

The  governing  equation  is  (3.5) 


(        +2\ 
\-r+ 

V  ) 


30       320       1    30        1     320 

+ +  ■ 


where 


&+      5S     r+dr*     Pel^ 


t0  -t  _+       r  u+  _  xl  rQ 


0  =  ^—,         r+=  — ,        x+ = 


t0-te  r0  RePr 


The  finite  difference  approximations  for  the  first  and  second  derivatives  with  error  of 
order  h  are  (figure  2.) 

90  0(,+U)  ~  0(,-u) 


3x+  (,j)  2hx 


+  o(hAx)  (4.1a) 


|*        =W^,-i)+Q(A2)  (41b) 


+  o(^2)  (4.1c) 
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d20  0itJ+l)  +  0{iJ_X)  -  20UJ)  | 


dr+~  (fj) 


h; 


+  <*ht) 


(4.  Id) 


When  these  approximations  are  plugged  into  the  governing  equation  (3.5),  the  equation 
becomes 


0U.j+D  +  %,j-\)  -  19(K. 


y«j) 


1 ft       a     ^ 


■*-*44a^uL)- 


(4.2) 


After  simplifying  the  expression  (4.2),  9{i })  becomes 


yM    i 


9n  „  =  ■ 


KPe2Dh2t      2hx      2hxJ 


yc-\.]) 


vPe2Dh2x      2hx      2hx 


'(/./+D 


1  1 

■  +  - 


K    2V«.,> 


i        i 


K   2V<,v 


)y 


(4.3) 


2  2 


2.  The  Boundary  Conditions 


The  boundary  conditions  are: 


a)       At  the  Center 


At  r  =0,  =  0 

dr+ 


(Symmetry  boundary  condition) 


1    d0 
When  this  boundary  condition  is  applied  to  the  governing  equation,  the  element 

r+  dr+ 

becomes  — .  After  applying  the  L'Hospital's  rule,  the  governing  equation  for  the  given 
boundary  condition  becomes 
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eo      cro_         d2e 


2    -,    +2 


8x+         dr+~      Pe~  dx+ 


(4.4) 


Since  the  temperature  distribution  at  the  center  on  the  radial  direction  is  symmetric, 
After  applying  the  finite  difference  method,  the  0(l  J}  at  the  boundary  becomes 


0. 


0+lj) 


1  1 


0 


KPe2Dh2x      2hXJ 


+  0, 


O-lj) 


f       1        |      1    ^ 
KPe2Dh2x+2hxJ 


+  0. 


(ij+l) 


(>J) 


+ 


p*dK    K 


b)       At  the  Wall 


(4.5) 


atr+=l. 


*<W)  =  ° 

for 

x+  >0 

80 

-0 

for 

x+  <0 

l5r+ 


(4.6) 


FOR  x+<0 


x+<0,       r+=l, 


80 
8r+ 


=  0 


The  governing  equation  for  the  given  boundary  condition  becomes 


820        1     8-0 


2     -s    +2 


8r+       PelD  8x+ 


=  0 


(4.7) 


80 
Since  the  wall  is  insulated,  — -  =  0 ,  0{i  .+1)  =  0{l  y  „  at  the  boundary. 

After  applying  the  finite  difference  approximations,  the  #(; ;)  at  the  boundary  becomes 
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e 


0"+W) 


°{i,j)  - 


<  1  ^ 

KP*lhlj 


+  0, 


('-!.;) 


f     1     ^ 
\Pelhl  j 


+  0 


(i,/-i) 


^  +  *,a 


(4.8) 


c)         As  X  -^   -  oo 


as  x+  ->-oo,         0(. ,.  — »  1 


(4.9) 


d)        Asx->  + 


00 


as  x+  -» +oo ,         0       ->  0 


(4.10) 


3.  The  Bulk  Mean  Temperature 


The  mean  temperature  at  a  position  "i"  is  (2.8) 


2      r°r 


(x)       Vr1   , 


hx,r)t(x,rfdr 


where  V  is  the  average  velocity  and  r0  is  the  radius  of  the  tube.  When    \t0u27wdr   is 

o 

subtracted  from  both  sides  and  dimensionless  parameter  r+  is  substituted  inside  the 
equation,  the  equation  becomes 

{tm-t0)=\(t-t0)u+2r+dr+ 


Let's  define  the  bulk  mean  temperature,   6m  -  y2 ^  > 

Ye  ~lo) 


(4.11a) 
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where    te    is    the    inlet    temperature    (as    x  — >  -oo )    and    to    is    the    constant    wall 
temperature.(Fluid  exit  temperature  as  x  — >  +00  ) 
Then  the  bulk  mean  temperature  9m  becomes 


6  -4ffl-r+2V        ,  r+dr+ 

m(x+)  IX  J    {x+,r+) 


(4.11b) 


6  is  found  out  from  the  numerical  calculation.  After  that  6m  is  calculated  by  solving  the 
equation  (4.1 1)  numerically. 


4.  The  Nusselt  Number 


The  Nusselt  number  is  (2.3) 


Nu  = 


hd 


When   the   heat   transfer   equations    q0=h(T0-Tm)    and    q0=-k 


dimensionless  parameter  r+  are  plugged  in  (2.3),  the  equation  becomes 


dr 


,    and   the 


r=R 


Nu  =  - 


2    30 


Omdr+ 


(4.12) 


rT=l 


The  finite  difference  approximation  for  the  equation  (4.12)  using  the  backward 
difference  expression  with  error  of  order  h    becomes 


Nu{l)  = 


e 


m(i) 


-  3>6(IJ)  +  4fl(,.ijH)  -  6{iJ_2) 

2/2.. 


(4.13) 


where  "j"  is  the  value  which  corresponds  to  r+=l . 
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5.  Heat  Absorption  Upto  the  Entrance 

Assume  a  differential  volume,   dV  =  Ijrrdr    between   x  =  -oo    and   x  =  0 .  The 
energy  stored  in  this  control  volume  is 

dE  =  cpAtdm  (4.14) 

where  At  is  the  temperature  difference  between  entrance  temperature  and  the  steady  state 
temperature  after  the  heating  starts,  Ar  =  t,xr\  -  te  and  dm   is  the  mass  of  the  control 

volume,  dm  =  pdV  .The  energy  absorbed  in  the  control  volume  becomes 

0  ro 

E absorbed  =  2^/* 'P  J   j(f  ~  te)rdrdx  .  (4.15) 

-00     0 

o  ro 
When      27zpcp  J  \t0rdrdx    is  subtracted  from  both  sides  and  (Vto)  is  divided  to  both 

-oo   0 

sides,  the  equation  (4.15)  becomes 

v  o  ro 

aWe*     =  In  J  \{G  -  \)rdrdx  (4.16) 

ih  ~to)PCp  -oo  0 

t  —  t  XI Y  Y 

where  6  = — .  As  the  non-dimensional  parameters  of  x+  =        °     and  y+  =  —  and 

their  first  derivatives  are  substituted  into  the  equation  (4.16),  the  equation  becomes 

F  °  ] 


'absorbed         _ 


.3 


(t0-te)PCpr0  -oo  0 


=  -27rPe  J  \(0  -  \)Y+dY+dx+  (4. 1 7) 
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B.    FLOW  INSIDE  THE  CIRCULAR  TUBE,  CONSTANT  HEAT 


FLUX  CASE 


1.  The  Governing  Equation 


The  governing  equation  is  (3.6) 


\  —  r 

v  J 


ee     dze     1  80      i    dle 

+ +• 


dx*      ^     r*dr*     Pel  ^ 


where 


0  =  ^^,        r+=^,         X+=^l± 
q0  d/k  r0  Red  Pr 


The  finite  difference  approximations  for  the  first  and  second  derivatives  with  the  error  of 
order  h2  are  (4.1a,b,c,d)  (figure  2.) 


d6_       =Gm,j)~G{i-\j) 
dx+  Pj)  2hx 


+  o(h;) 


_*ft/*i)      e«J-±  +  0{h2y) 


dr*  oj) 


2h 


d20  0V*U)+0V-u)-20<tj) 


dx+   (ij 


u) 


hi 


+  o(h;) 


dr+  (ij)  hy 
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When  these  approximations  are  plugged  into  the  governing  equation  (3.6),  the  equation 
becomes 


1   (ec*ij)  +  9o-ij)-20oj) 


Pelr 


6U.j+\)  +  e(:j-D  ~  29U.J) 


■Vcj) 


0(:j+D     9(>.j-\) 


2hy 


4-yh{^^)=o 


(4.18) 


After  simplifying  the  expression  (4.18),  6{iJ)  becomes 


o„  n  = 


L^+yLi_± 


^Pe2Dh2x      2hx      2hx 


+  *<.-.. 


(i-U) 


i     y,.n .  i 


KPe2Dhl      2hx      2hx 


+e, 


C.;+i) 


I_  1 

vi  2h>y«.»j 


( \       i 


'o.j-d 


Khl      2VW) 


(4.19) 


2  2 

^r-^r  +  ^r 


2.  The  Boundary  Conditions 

The  boundary  conditions  are: 


a)        At  the  Center 


At  r+=0. 


30 


=  0 


(Symmetry  boundary  condition) 


1    30 

When  this  boundary  condition  is  applied  to  the  governing  equation,  the  element 

r+  dr+ 

becomes  — .  After  applying  the  L'Hospital's  rule,  the  governing  equation  for  the  given 
boundary  condition  becomes 


d6     „  d20        1     d26 

=  2 r  + 


dx+         5r+*      Pe-  dx+ 


(4.20) 
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Since  the  temperature  distribution  at  the  center  on  the  radial  direction  is  symmetric, 
After  applying  the  finite  difference  method,  the  6(l  y)  at  the  boundary  becomes 


0 


0+U) 


0, 


(i.7) 


Pe2 h2      2h 


+  &, 


0-1,;) 


1  1 

+ 


\Pehhx 


+  6 


(',;+D 


Pe2Dh2x     h; 


b)        At  the  Wall 


(4.21) 


atr+=l. 
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0 


1 

for 

x+  >0 

z 

for 

x+  <0 

(4.22) 


FORx+>0i 


r    = 


•    u  a#      i 

so  the  equation  becomes =  — . 

dr+         2 


The  finite  difference  approximation  for  this  equation  is 


fyfj+i)     *(fj-i)         1 


2/2. 


When  the  boundary  condition  is  applied  to  the  governing  equation  (3.6),  the  equation 
becomes 


d20        1  1     d20 


+  ■ 


2    ~    S 


dr+       2r+      Pe1  dx+ 


=  0 


(4.23) 
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As  the  finite  difference  approximations  are  substituted  in  the  equation  above,  the  equation 
becomes 


ft 


('■+!,;) 


0, 


<       1       ^ 
KPeWxJ 


+  0. 


(H;) 


f      1       ^ 

v^A2y 


+  0, 


(u-i) 


(»■,/) 


f^ 


v  y  J 


J_    J. 

^     2 


PejAj  +  A* 


(4.24) 


FOR  x  <0 


x+<0,       r+=l, 


30 
dr+ 


=  0 


The  governing  equation  for  the  given  boundary  condition  becomes 


d20         1      d20      A 

+ : -  =  0 


dr+1      PeD  dx 


,2    .    ,2 


(4.25) 


Since  the  wall  is  insulated,  — -  =  0 ,  0{i  j+l)  =  6{l  ;_1}  at  the  boundary. 

After  applying  the  finite  difference  approximations  to  (4.25),  the  6(i })  at  the  boundary 
becomes 


^(<+w; 


0OJ)  - 


f      i      N\ 


KPelhlj 


r     ,      \ 


+  0. 


(i-lj) 


Pe2Dh2xJ 


+  0, 


0,7-1) 


Kh'h 


+ 


Pe2Dh;     h; 


(4.26) 


c)        As  x  -> 


CO 


as  x+  ->-oo,         0{i  n  ->0 


(4.27) 


30 


d)        As  x  ->  + 


00 


as  x    -»  +oo 


°(U)  ">  "2*+  + 


4  ,2^1 


7 


+ 


48       8         2 


v 


(4.28) 


3.  The  Bulk  Mean  Temperature 

The  mean  temperature  at  a  position  "i"  is  (2.8) 


2      r°r 


f™=V2     \UtVdr 

Vr0    o 


where  V  is  the  average  velocity  and  r0  is  the  radius  of  the  tube.  When    \teu2izrdr  is 

o 

subtracted  from  both  sides  and  dimensionless  parameter  r+  is  substituted  inside  the 

equation,  the  equation  becomes 


{tm-te)=\(t-te)u+2r+dr+ 
o 

Let's  define  the  bulk  mean  temperature, 

fc-0 


ft. 


q'od/k 


(4.29) 


(4.30a) 


where  te  is  the  inlet  temperature  (as  x  — >  -co ). 
Then  the  bulk  mean  temperature  6m  becomes 


ft 


•4-sy 


dr 


(4.30b) 


31 


6  is  found  out  from  the  numerical  calculation.  After  that  6m  is  calculated  by  solving  the 
equation  above  numerically. 

4.  The  Nusselt  Number 

The  Nusselt  number:  Nu  =  — . 

k 

The  heat  transfer  from  the  wall: a'  =  hit,  t  +  ,  -  O  • 

10  v  (x  ,r  =1)  m ' 

The  temperature  at  the  position  of  (x+,r+=l):  t  ,  +      =te-6q"0  d/k 

The  average  temperature  at  (x+):  tm  -te-  0mqo  d/k 

When  these  equations  are  substituted  into  the  Nusselt  number  equation,  Nusselt 
number  becomes: 

Nu  =  — l- —  (4.31) 

0M-O 

ffl 

5.  Heat  Absorption  Upto  the  Entrance 

Assume  a  differential  volume,   dV  =  Irurdr   between  x  =  — oo   and  x  =  0 .  The 
energy  stored  in  this  control  volume  is  (4. 14) 

dE  =  cpAtdm 

where  At  is  the  temperature  difference  between  entrance  temperature  and  the  steady  state 
temperature  after  the  heating  starts,  At  =  f(x  r)  -  te  and  dm   is  the  mass  of  the  control 

volume,  dm  =  pdV  .The  energy  absorbed  in  the  control  volume  becomes 
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0   ro 

E absorbed  =  ^Pc  P  j   J(* "  te)rdrdx  .  (4.32) 


p 

-oo   0 


When  both  sides  are  divided  by  q0  d/k  ,  the  equation  (4.32)  becomes 


F  o  ro 

absorbed       =  _A?J.   J     J^^  (4  33) 


^loPCpro/k  -»  0 

f  —  /  x/r  a* 

where  #  =  — £ .  As  the  non-dimensional  parameters  of  x+  =        °    and  r+  =  —  and 

q0d/k  RePr  r0 

their  first  derivatives  are  substituted  into  the  equation  (4.33),  the  equation  becomes 

E absorbed       =  _Anpe  J  J#^+aX  +  (4-34) 


q"oPcpr?/k 


C.    FLOW  BETWEEN  THE  PARALLEL  PLATES,  CONSTANT 
WALL  TEMPERATURE  CASE 
1.  The  Governing  Equation 

The  governing  equation  is  defined  in  equation  (3.1 1)  as 

1      d26      A  d26      u+  d6 

+  4- 


*k  dx*2        dy+2       2  &+ 


where 


t0-te       '        h  RePr  h 

The  finite  difference  approximations  for  the  first  and  second  derivatives  with  error  of 
order  h  are  (4.1a,b,c,d)  (figure  2.) 
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dx+  (,,;)  2hx 


o(h;) 


dr+{1J)  2/2 


dx+2  «,„  h2x 


o(h;) 


d2&         _e(ij+i)+0(i,j-i)-20V,fi 


dr+~  (,j) 


h; 


+  o(h2v) 


When  these  approximations  and  u+  are  plugged  into  the  governing  equation  (3.1 1),  6{l  n 


becomes 


'o+i.j) 


i      fyij)     3 


KPe2Dh2x        8hx       &hx 


O-IJ) 


1  3-^.;>.  +  .  3 


Pel.hl       ihx       %hx 


+  0, 


(>,;+D 


f      \ 


Sj 


r  .  \ 


'{>j-\) 


Khyj 


(4.35) 


'(•j) 


8 


Pe\hl      K 


2.  The  Boundary  Conditions 


The  boundary  conditions  are: 


a)        At  the  Center 


At  y+=0,      j=l, 


80 


=  0 


(Symmetry  boundary  condition) 
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When  this  boundary  condition  is  applied  to  the  governing  equation,  the  governing 
equation  becomes 


1      d26      A  d26      3  80 
— +  4- 


Pei  dx+'        dy+~      4  dx+ 


(4.36) 


Since  the  temperature  distribution  at  the  center  on  the  radial  direction  is  symmetric, 


^(«.y+i)  ~  *qj-i) 


at  the  center. 


After  applying  the  finite  difference  method  to  (4.36),  the  6{i  n  at  the  boundary  becomes 


e, 


(i+i.D 


dU.D  = 


yPe~Dh;     %hxJ 


+  0, 


('-i,i) 


1  3 

+ 


yPe2Dh;     8/^ 


+  0, 


(',2) 


+ 


PeoK     K 


(4.37) 


b)        At  the  Wall 


aty+=l, 


dO 


dy+ 


=  0 


for      x+  >  0 


for      x+  <  0 


(4.38) 


FOR  x+<0 


x+<0,       y+=l, 


dy+ 


=  0 


The  governing  equation  for  the  given  boundary  condition  becomes 

1      d29     A  d26      n 

: r  +  4 r  =  0 

Pe~Dh  dx+'        dy+~ 


(4.39) 
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d0 
Since  the  wall  is  insulated,  =  0,  6{i  +1)  =  6{l  _1}  at  the  boundary. 

dy+ 

After  applying  the  finite  difference  approximations,  the  0{i ;)  at  the  boundary  becomes 


^(<+U) 


°{u)  - 


KPelhl  J 


+  0 


(i-hj) 


f        1        ^ 

\p*2Dhl  j 


+  *<U-i] 


2  8 

+  ■ 


^2   K 


(4.40) 


e)        As  x  -> 


00 


as  x+^-oo,         #(1;)  — »1 


(4.41) 


f)         As  x  -»  + 


CO 


as  x+  ->+oo,         #(lJ)  -»0 


(4.42) 


3.  The  Bulk  Mean  Temperature 

The  mean  temperature  is  defined  in  equation  (2.9)  as 

When  the  dimensionless  parameters  u+,  y+  are  substituted  in  (2.9),  the  tm  becomes 


tm  =  \u+tdy+ 
o 

where  V  is  the  average  velocity  and  h  is  the  distance  from  the  centerline  to  the  plate. 

1 
When  t0  =  jt0u+dy+  is  subtracted  from  both  sides,  the  equation  becomes 

0' 
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(to-tm)=)(t0-t)u+dy+ 

o 

Let's  define  the  bulk  mean  temperature, 

*"=fef  (4-43a) 

where    U    is    the    inlet    temperature    (as    x  — >  -co )    and    to    is    the    constant    wall 
temperature. (Fluid  exit  temperature  as  x  — >  +00 ) 
Then  the  bulk  mean  temperature  0m  becomes 


em=^){l-y+2^y\  (4.43b) 


6  is  taken  from  the  numerical  calculation.  After  that  0m  is  calculated  by  solving  the 
equation  (4.43)  numerically. 

4.  The  Nusselt  Number 

The  Nusselt  number  is  defined  as  (2.3) 

Nu  =  — h-. 

k 

30 
When  the  heat  transfer  equations  q0  =  h(T0  -Tm),  qQ  =  k — 

dy 


,  and  the  dimensionless 

y=h 


parameter  r+  are  plugged  in  (2.3),  the  equation  becomes 

4    60 


Nu  =  - 


0mdy+ 


(4.44) 
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The  finite  difference  approximation  for  the  equation  (4.44)  using  the  backward  difference 
expression  with  error  of  order  h    becomes 


Nu,  = 


f 


0mi 


2hy 


(4.45) 


5.  Heat  Absorption  Upto  the  Entrance 

Assume  a  differential  volume,   dV  =  bdxdy   between  x  =  -oo   and  x  -  0 .  The 
energy  stored  in  this  control  volume  is  (4.14) 

dE  =  cpAtdm 

where  At  is  the  temperature  difference  between  entrance  temperature  and  the  steady  state 
temperature  after  the  heating  starts,  At  =  t^x  ^  -  te  and  dm   is  the  mass  of  the  control 

volume,  dm  =  pdV  .The  energy  absorbed  in  the  control  volume  becomes 

^absorbed  =  *f*p    1    J('  "  O^  •  (4-46) 

-oo-A 


0  h 

'P 


When  2bpcp  J  \t0dydx  is  subtracted  from  both  sides  and  (te-t0)  is  divided  to  both 


-oo  0 


sides,  the  equation  (4.46)  becomes 


F  o  h 

=  2  J  \{6-\)dydx  .     (4.47) 


bpc(te-t0)       ^o 
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where     0  = — .    As    the    non-dimensional    parameters    of    x+  =  — - — -     and 

r+  =  —  and  their  first  derivatives  are  substituted  into  the  equation  (4.47),  the  equation 


becomes 


F  °  1 


'abs°'bed       =  -APe  j  \(6  -  l)dy+dx+  (4.48) 


hbpc(t0-te)  .ooo 


D.  FLOW  BETWEEN  THE  PARALLEL  PLATES,  CONSTANT 
HEAT  FLUX  CASE 

1.  The  Governing  Equation 

The  governing  equation  is  (3.12) 

1      d26      A  d26      u+  d6 

+  4- 


Pe%,  dx*2        dy+2       2  &+ 


where 


■0  =  4^,      y*^,       x*=^^,       D„=4h 
q.DJk'       y       h'  RePr 

The  finite  difference  approximations  for  the  first  and  second  derivatives  are  (4.1a,b,c,d) 

dx+  (,,;)  2hx 


d0_       =  A(W     °(u-n+(Kh2) 

dr+{lJ)  2hy 


39 


d2G  0(,+u)+^-u)-2*o,,) 


dx+  a  j) 


+  o(h2) 


d20 

dr+2{l 


0w+\"-W™  +  0(h2) 


i'j) 


When  these  approximations  and  u+  are  plugged  into  the  governing  equation,  the  bulk 
temperature  at  the  position  (i  J)  is  written  as 


e, 


IMJ) 


^)=" 


1      J^       3  ) 


^Pe2Dh]       %hx       %htJ 


+8, 


(i-U) 


i      H, 


+  0, 


('.y+D 


+e, 


Oj-U 


2  8 

■  +  - 


Pelh2x     h; 


(4.49) 


2.  The  Boundary  Conditions 


The  boundary  conditions  are: 


a)        At  the  Center 


At  y+=0,      j=ls 


de 

dy+ 


=  0 


(Symmetry  boundary  condition) 


When  this  boundary  condition  is  applied  to  the  governing  equation  (3.12),  the  governing 
equation  becomes 


i    d2e    A  d2e    3  de 

+  4- 


Pe-Dh  dx+         dy+~      4  dx+ 


(4.50) 


Since  the  temperature  distribution  at  the  center  on  the  radial  direction  is  symmetric, 
0<*j+D  =evj-r)  at  the  center. 
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After  applying  the  finite  difference  method,  the  9{i ;)  at  the  boundary  becomes 


e, 


(7+1,1) 


0(i,D  ~ 


KPe2Dh2x      ShxJ 


+  9, 


('-i,i) 


f      1       |     3  % 
KPe2Dh2x  +  Mx  j 
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(»■,;+!) 


^4^2  +  ^ 


(4.51) 


b)        At  the  Wall 


aty=l,      = 

dy+ 


0 
J_ 

4 


/or      x+  <  0 


for      x+  >  0 


FOR  x+<0 


x+<0,         y+=l,  ^.  =  0 

6>+ 


The  governing  equation  for  the  given  boundary  condition  becomes 

i    d2e    .  d2e    . 

Pel  dx+         dy+ 


(4.52) 


Since  the  wall  is  insulated,  — -  =  0 ,  9{l  J+l)  =  9{i ,_»  at  the  boundary. 


After  applying  the  finite  difference  approximations,  the  9{i  n  at  the  boundary  becomes 


6, 


0+1,./) 


e^j)  ~ 
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KPelhlj 


+  *<m. 
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KP'lhlj 


0(U-D 
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Pe2Dh2x  +  h) 


(4.53) 
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FOR  x+  >  0 


x+>0,      y+=l, 


60         J. 

dy+  '      4 


The  governing  equation  for  the  given  boundary  condition  becomes 

i    d2e    A  d2e    n 

+  4 -  =  0 


Pe-Dh  dx+         dy+ 


(4.54) 


d0_     _]_        Qqj+i)  ~ &(u-i) 
dy+  '      4 '  2/z 


1  h 

- ,         0{lJ+X)  =  6{lJ_x)  -—      at  the  boundary. 


After  applying  the  finite  difference  approximations,  the  6{l ;)  at  the  boundary  becomes 


e. 


('"+1.7) 


eiKi)   ~ 


J_ 

KPelhlj 


+  0, 


0-1.7) 


_1_ 

yPe2Dh2xJ 


e, 
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Pe2Dh;     h; 


(4.55) 


c)        As  x  -> 


oo 


as  x+  ->  -oo ,         6U  n  ->  0 


(4.56) 


d)        As  x  -»  + 


00 


as  x    ->+oo 


*(/j>  =  ~2x+  + 


39        3 


(    .4 


1120     16 


y        +2 
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(4.57) 


3.  The  Bulk  Mean  Temperature 

The  mean  temperature  is  (2.9) 
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When  the  dimensionless  parameters  u+,  y+  are  substituted  in  (2.9),  the  tm  becomes 


tm  =  )u+Tdy+ 
o 

where  V  is  the  average  velocity  and  h  is  the  distance  from  the  centerline  to  the  plate. 

1 
When  te  =  \teu+dy+  is  subtracted  from  both  sides,  the  equation  becomes 

o 

{tm-te)=){t-te)u+dy+ 
o 

Let's  define  the  bulk  mean  temperature, 

^=%T77>  (4.58a) 

where  te  is  the  inlet  temperature  (as  x  — >  -oo  ).  Then  the  bulk  mean  temperature  0m 
becomes 

^  0 

6  is  taken  from  the  numerical  calculation.  After  that  6m  is  calculated  by  solving  the 
equation  (4.58)  numerically. 


4.  The  Nusselt  Number 


hD 
The  Nusselt  number:  Nu  =  — - 

k 


The  heat  transfer  from  the  wall: q0  =  h(t  +  ^     -tm) 


The  temperature  at  the  position  of  (x+,r+=l):  t  +  r+_    =te  -  0qo  Dh/k 
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The  average  temperature  at  (x+):  tm  -te-  6mq0  Dh/k 

When  these  equations  are  substituted  into  the  Nusselt  number  equation,  Nusselt  number 
is  calculated  as 

M»--i-  (4.59) 

m 

5.  Heat  Absorption  Upto  the  Entrance 

Assume  a  differential  volume,   dV  =  bdxdy  between  x  =  -oo   and   x  =  0 .  The 
energy  stored  in  this  control  volume  is  (4.14) 

dE  =  cpAtdm 

where  At  is  the  temperature  difference  between  entrance  temperature  and  the  steady  state 
temperature  after  the  heating  starts,  At  =  t^x  y)  -  te  and  dm  is  the  mass  of  the  control 

volume,  dm  =  pdV  .The  energy  absorbed  in  the  control  volume  becomes 

^absorbed  =  *>fXp    1  +J('  ~  Qbffydx  .  (4.60) 

-oo  — h 

When  both  sides  are  divided  by  q0  Dh/k ,  the  equation  (4.60)  becomes 

F  oh 

pbssrbed      =  _g  J    j$dydx  (4.61) 

VobpCp  hlk         -«  o 

t  —t  Ixl  D  r 

where  6  =  — ^ .  As  the  non-dimensional  parameters  of  x+  =  — - — -  and  r+  =  —  and 

q0Dhlk  RePr 

their  first  derivatives  are  substituted  into  the  equation  (4.61),  the  equation  becomes 

F  °  1 

'absorbed   =  _\6pe   j   fay* dx+  (4.62) 


q0bh  ja  _k0 
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E.  THE  EFFECTS  OF  WALL  CONDUCTION,  FLOW  INSIDE  THE 
CIRCULAR  TUBE,  CONSTANT  WALL  TEMPERATURE  CASE 
1.  The  Governing  Equation 

The  governing  equation  is  (3.5) 


where 


e  =  ^L 


h-K 


r    =  —  ,        x    = 


x/r0 


RePr 


The  finite  difference  approximations  for  the  first  and  second  derivatives  with  error  of 
order  h  are  (4.1  a,b,c,d)  (figure  2.) 


de_     J(MJ)   °(hj)+0(^) 

dx+  (,,;■)  2hx 


+  o(h2y) 


d20        =9v+ij)+0(>-ij)-20(u) 
dx+2  iu)  h2x 


+  o(h;) 
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When  these  approximations  ai    plugged  into  the  governing  equation  (3.5),  the  equation 
becomes 


1      (d0*lj)+d0-lj)~2d(:.j) 


Pel 


@(i.j+V  +  &0,j-\)  ~  2&U.J) 


yoj) 


\ij+D    fyi.j-n 


2h 


y         ) 


\-yl\ 


2       1   %+\.j)  ~  &(i-l.j) 


2k. 


(4.63) 


After  simplifying  the  expression  (4.63),  the  equation  becomes 


*0,>  = 


i    ,  y,j)    i 

Pelti      2h.      2h 


1       y,n  .    1 


<'-'•-"      D^,2 


+  *(..,♦.) 


1  1 

— r  +  - 


'l  1 


v^;    2Vo,)y 


'c.;-D 


^;    2V..,>; 


(4.64) 


2  2 


2.  The  Boundary  Conditions 


The  boundary  conditions  are: 


a)        At  the  Center 


At  r+=0, 


80 

dr  + 


=  0 


(Symmetry  boundary  condition) 


1    dO 
When  this  boundary  condition  is  applied  to  the  governing  equation,  the  element 

r+  dr+ 

becomes  — .  After  applying  the  L'Hospital's  rule,  the  governing  equation  for  the  given 


boundary  condition  becomes 


36      „  d26        1     d26 

=  2 -  + 


dx+         dr+'      Pe~  dx+' 


(4.65) 


Since  the  temperature  distribution  at  the  center  on  the  radial  direction  is  symmetric, 
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After  applying  the  finite  difference  method,  the  0(!  ,»  at  the  boundary  becomes 


0, 


O+lj) 


0 


iUJ) 


Pp2 h2      2h 

\reDnx        ^nx  J 


+  0, 


O-iJ) 


i       i 

Pp2  h2       2h 
\reDnx        *-nx  j 


+  0, 


(',;+!) 


r4y 


Pe2Dh2  +  h] 


b)        At  the  Wall 


(4.66) 


atr+=l, 


@wal!  -  &  fluid 
dewall   _  de fluid 


{  dr+  dr+ 


The  heat  fluxes  through  the  fluid  and  the  wall  are  equal  to  each  other  at  r+=l. 


K  ^ solid    _   7.  ^ fluid 

*  solid     ^      ~K flwd     ^+ 


>nd 


Using  the  2      order  forward  and  backward  difference  method,  the  equality 
becomes 


30{iJ)  -40(ij-i)  +0(,j-2)  I     ,  (-w(u)  +  4%y+i)  -0d,j+ 


2/2 


yf 


=  k. 


1-2) 


2/z 


ys 


(4.67) 


If  hys=hyf>  the  bulk  temperature  at  the  fluid- wall  boundary  becomes 


0, 


%j+l)  +  k  fyj-1) 


kf 

e(ij+2)+-reu 

Ks 


\Uj-2) 


(W) 


(4.68) 
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c)        On  the  Outer  Edge  of  the  Wall 


FOR  x+<0 


For  x+<0  and   r+  =  1  +  dy ,     the  outer  surface  of  the  solid  wall  is  insulated  where 


Sy  = 


solid  _  wall  _  thickness 
tube    inner    radius 


The  governing  equation  for  the  given  boundary  condition  becomes 

d2e      1   d2e 


-r  + 


2    ^    +2 


dr+'      Pe~  dx 


=  0 


(4.69) 


dO 
Since  the  wall  is  insulated,  — -  =  0 ,  9{l  J+l)  =  0(i  ;._15  at  the  boundary. 

After  applying  the  finite  difference  approximations,  the  6{l  n  at  the  boundary  becomes 


9 


('+1,7) 


°{ij)  ~ 


\Pehh*  J 


+  9 


(«■-!.;) 


Pe2Dh2xJ 


+  0, 


0,;-i) 


+ 


Pe2Dh2x      h2y 


(4.70) 


FORx+>l 


°iu)  =  0 


(4.71) 


d)        As  x  ->  -  oo 


as  x+  -»  -co ,         0(i  .j  — »  1 


(4.72) 


e)        As  x  -»  + 


CO 
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asx+->+oo,  0(,-.;)-»°  (4-73) 


3.  The  Bulk  Mean  Temperature 

The  mean  temperature  at  a  position  "i"  is  (2.8) 


2    r° 

VrQ      0 


Tr 


where  V  is  the  average  velocity  and  r0  is  the  radius  of  the  tube.  When    \t0ul7irdr  is 

o 

subtracted  from  both  sides  and  dimensionless  parameter  r+  is  substituted  inside  the 

equation,  the  equation  becomes 

{tm-t0)=\{t-t0)u+2r+dr+ 
o 

Let's  define  the  bulk  mean  temperature,  0m  -  -V2 2/,  where  te  is  the  inlet  temperature 

V*  e  —  *  o  ) 

(as  x  -»  -oo )  and  to  is  the  constant  wall  temperature.(Fluid  exit  temperature  as  x  ->  +oo ) 
Then  the  bulk  mean  temperature  Om  becomes 


^-w=4jt-r+aW+  <4-74) 


^  is  found  out  from  the  numerical  calculation.  After  that  6m  is  calculated  by  solving  the 
equation  (4.74)  numerically. 
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4.  The  Nusselt  Number 

The  Nusselt  number  is  defined  as  (2.3) 

Nu  =  — . 
k 

When    the    heat    transfer    equations    q0  =  h(t0  -tm)    and    q0  =  -k 
dimensionless  parameter  r+  are  plugged  in  (2.3),  the  equation  becomes 


30 

dr 


,    and    the 


r=R 


Nu  =  - 


2    d6 


emdr+ 


r+=l 


(4.75) 


The  finite  difference  approximation  for  the  equation  (4.75)  using  the  backward 
difference  expression  with  error  of  order  h    becomes 

2    f-S^,^,,,,.,,-^,, 


Nu(n  = 
(0      e 


m(i) 


2/2, 


(4.76) 


where  "j"  is  the  value  which  corresponds  to  r+=l, 


5.  Heat  Absorption  Upto  the  Entrance 

Assume  a  differential  volume,   dV  =  27irdr    between  x  =  -co   and  x  =  0 .  The 
energy  stored  in  this  control  volume  is  (4. 14) 

dE  =  cpAtdm 

where  At  is  the  temperature  difference  between  entrance  temperature  and  the  steady  state 
temperature  after  the  heating  starts,  At  =  t^x  r)  -  te  and  dm   is  the  mass  of  the  control 

volume,  dm  =  pdV  .The  energy  absorbed  in  the  control  volume  becomes 
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o  ro 
E absorbed  =  2^P°P  J   }('  ~  te)rdrdx  .  (4.77) 

-oo    0 
0    ro 

When      Inpc    J    \t0rdrdx    is  subtracted  from  both  sides  and  (Vto)  is  divided  to  both 

-oo   0 

sides,  the  equation  (4.77)  becomes 

p  o  r0 

— absorbed     =  In  \  \{9  -  X)rdrdx  (4.78) 

(te  ~*o)PCp  -oo   0 

t  —  t  XI Y  Y 

where  6  = — .  As  the  non-dimensional  parameters  of  x+  =        °     and  y+  =  —  and 

te~t0  RePr  y0 

their  first  derivatives  are  substituted  into  the  equation  (4.16),  the  equation  becomes 

F  o  i 

absorbed      =  -2nPe  \  \{6  -  X)Y+dY+dx+  (4.79) 

(t0-te)pcr0  _oo  o 
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V.       NUMERICAL  RESULTS  AND  DISCUSSION 


A.    BULK  MEAN  TEMPERATURE  VS.  X 

The  bulk  mean  temperature  graphs  of  the  fully  developed  laminar  flow  inside  the 
circular  tube  and  parallel  plates  for  constant  wall  temperature  and  heat  flux  cases  were 
plotted  in  figures  8-11.  The  Peclet  numbers  0.5,  1,  2,  3,  5,  10,  50,  and  100  were  used  for 
the  plots  to  see  how  the  Peclet  number,  Pe=RePr,  affects  the  axial  conduction  in  the  flow. 
As  it  is  seen  from  the  figures,  when  the  Peclet  number  gets  smaller,  the  axial  conduction 
effects  of  the  fluid  becomes  very  important  and  the  temperature  of  the  fluid  upto  the 
entrance  increases.  The  axial  conduction  of  the  fluid  has  to  be  considered  for  these  cases 
and  must  not  be  neglected  for  the  Peclet  numbers  less  than  about  10.  When  the  Peclet 
number  gets  larger,  the  effects  of  the  axial  conduction  decrease.  As  it  is  seen  from  the 
figures  8-11,  for  the  Peclet  numbers  50  and  100,  the  bulk  mean  temperature  of  the  fluid  at 
the  position  x+=0  is  almost  equal  to  the  entrance  temperature,  which  means  that  there  is 
almost  no  heat  absorption  or  temperature  increase  upto  the  entrance.  The  axial  conduction 
of  the  fluids  can  be  neglected  for  these  high  Peclet  numbers. 

1.    Flow  Inside  the  Circular  Tube,  Constant  Wall  Temperature 

Case 

The  graphs  of  the  bulk  mean  temperature  vs.  x  for  the  circular  tube,  constant  wall 
temperature  flow  was  plotted  for  various  Peclet  numbers  in  figure  8. a  and  8.b  where  theta 
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bulk  mean  is    Om  =  V5 — ^l  and  the  nondimensional  axial  distance  x+  is  x+  = 


fc-0 


RePr 


The  validity  of  the  figure  8.b  has  been  checked  by  comparison  with  Ref.5. 


CIRCULAR  TUBE  CONSTANT  WALL  TEMP.  FLOW.  THETA  BULK  MEAN 


Figure  8.a  The  axial  bulk  mean  temperature  distribution  for  circular  tube  constant  wall  temperature  flow 


CIRCULAR  TUBE  CONSTANT  WALL  TEMP.  FLOW.  THETA  BULK  MEAN 

■*e- 


Figure  8.b  The  enlarged  view  of  figure  8. a 
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2.    Flow  Inside  the  Circular  Tube,  Constant  Heat  Flux  Case 


The  graphs  of  the  bulk  mean  temperature  vs.  x  for  the  circular  tube,  constant  heat 
flux  flow  was  plotted  for  various  Peclet  numbers  in  figure  9. a  and  9.b  where  theta  bulk 

mean  is  0„  =  M; — —  and  the  nondimensional  axial  distance  x+  is  x+  =  ■ 


q0dlk 


Re^Pr 


CIRCULAR  TUBE  CONSTANT  HEAT  FLUX  FLOW.  THETA  BULK  MEAN 
0, ,0 rA- 


Figure  9.a  The  axial  bulk  mean  temperature  distribution  for  circular  tube  constant  heat  flux  flow 


CIRCULAR  TUBE  CONSTANT  HEAT  FLUX  FLOW,  THETA  BULK  MEAN 


-0.5       -0.4       -0.3       -0.2       -0.1  0  0.1         0.2        0.3         0.4         0.5 


Figure  9.b  The  enlarged  view  of  figure  9. a 
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3.    Flow  Between  the  Parallel  Plates,  Constant  Wall  Temperature 
Case 

The  graphs  of  the  bulk  mean  temperature  vs.  x+  for  the  parallel  plates,  constant  wall 

temperature  flow  was  plotted  for  various  Peclet  numbers  in  figure  lO.a  and  lO.b  where 
theta  bulk  mean  is  q   -  V°  ~tm>  and  the  nondimensional  axial  distance  x+  is 


fc-0 


x    = 


RePr 


PARALLEL  PLATES  CONSTANT  WALL  TEMP    FLOW,  THETA  BULK  MEAN 


Figure  lO.a  The  axial  bulk  mean  temperature  distribution  for  parallel  plates  constant  wall  temperature  flow 


PARALLEL  PLATES  CONSTANT  WALL  TEMP    FLOW.  THETA  BULK  MEAN 
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Figure  lO.b  The  enlarged  view  of  figure  10. a 
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4.    Flow  Between  the  Parallel  Plates,  Constant  Heat  Flux  Case 


The  graphs  of  the  bulk  mean  temperature  vs.  x  for  the  parallel  plates,  constant  heat 
flux  flow  was  plotted  for  various  Peclet  numbers  in  figure  1 1  .a  and  1 1  .b  where  theta  bulk 

mean  is  6„  =  ^4 —  and  the  nondimensional  axial  distance  x+  is  x+  = . 


q0Dh/k 


RePr 


PARALLEL  PLATES  CONSTANT  HEAT  FLUX  FLOW.  THETA  BULK  MEAN 
0, -e rA- 


Figure  ll.a  The  axial  bulk  mean  temperature  distribution  for  parallel  plates  constant  heat  flux  flow 


PARALLEL  PLATES  CONSTANT  HEAT  FLUX  FLOW.  THETA  BULK  MEAN 

0, 1 » *- 
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Figure  ll.b  The  enlarged  view  of  figure  1  l.a 
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B.    NUSSELT  NUMBER  VS.  X 

The  nusselt  number  graph  of  the  fully  developed  laminar  flow  inside  the  circular 
tube  and  parallel  plates  for  constant  wall  temperature  and  heat  flux  cases  were  plotted  in 
figures  12-15.  The  Peclet  numbers  0.5,  1,  2,  3,  5,  10,  50,  and  100  were  used  for  the  plots 
to  see  how  the  Peclet  number,  Pe=RePr,  affects  the  nusselt  number  distribution  in  the 


flow. 


1.    Flow  Inside  the  Circular  Tube,  Constant  Wall  Temperature 


Case 


The  graph  of  the  nusselt  number  vs.  x  for  the  circular  tube,  constant  wall 
temperature  flow  was  plotted  for  various  Peclet  numbers  in  figure  12  where  x+  is 


x    = 


2xld 
RePr 


CIRCULAR  TUBE  CONSTANT  WALL  TEMP.  FLOW.  NUSSELT  NUMBER 


Figure  12.  The  nusselt  number  distribution  for  circular  tube  constant  wall  temperature  flow 
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The  nusselt  numbers  approach  to  the  values  in  those  of  table  1  when  x  gets  larger  as  seen 
in  figure  12. 


Pe 

0.5 

1 

2 

3 

5 

10 

50 

100 

Nu.  no 

4.0971 

4.0276 

3.9224 

3.8506 

3.767 

3.6948 

3.6586 

3.6572 

Table  1.  The  nusselt  number  for  various  Peclet  numbers  for  circular  rube  constant  wall  temperature  case  as 
x  goes  to  infinity 


2.    Flow  Inside  the  Circular  Tube,  Constant  Heat  Flux  Case 

The  graph  of  the  nusselt  number  vs.  x  for  the  circular  tube,  constant  heat  flux  flow 


was  plotted  for  various  Peclet  numbers  in  figure  13  where  x+  is  x+  = 


Re^Pr 


CIRCULAR  TUBE  CONSTANT  WALL  TEMP.  FLOW,  NUSSELT  NUMBER 


Figure  13.  The  nusselt  number  distribution  for  circular  tube  constant  heat  flux  flow 


The  nusselt  number  is  independent  from  the  Peclet  number  for  this  case,  and  approaches 


to  4.364  where  x  goes  to  infinity. 


59 


3.    Flow  Between  the  Parallel  Plates,  Constant  Wall  Temperature 


Case 


The  graph  of  the  nusselt  number  vs.  x  for  the  parallel  plates,  constant  wall 
temperature  flow  was  plotted  for  various  Peclet  numbers  in  figure  14  where  x+  is 

.      2x1  Dh    ^  .    „ 

x+  = - ,  Dh  is  4h. 

RePr 


PARALLEL  PLATES  CONSTANT  WALL  TEMP.  FLOW.  NUSSELT  NUMBER 
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Figure  14.  The  nusselt  number  distribution  for  parallel  plates  constant  wall  temperature  flow 


The  nusselt  numbers  approach  to  the  values  in  those  of  table  2  when  x  gets  larger  as  seen 


in  figure  12. 


Pe 

0.5 

1 

2 

3 

5 

10 

50 

100 

Nu.  no 

8.0588 

8.0059 

7.9168 

7.8464 

7.7471 

7.6303 

7.5457 

7.542 

Table  2.  The  nusselt  number  for  various  Peclet  numbers  for  parallel  plates  constant  wall  temperature  case 
as  x  goes  to  infinity 


60 


3.    Flow  Between  the  Parallel  Plates,  Constant  Heat  Flux  Case 

The  graph  of  the  nusselt  number  vs.  x  for  the  parallel  plates,  constant  heat  flux 

flow  was  plotted  for  various  Peclet  numbers  in  figure  15  where  x+  is  x+  =  -      -1- ,  Dh  is 


4h. 


RePr 


PARALLEL  PLATES  CONSTANT  HEAT  FLUX  FLOW,  NUSSELT  NUMBER 


Figure  15.  The  nusselt  number  distribution  for  parallel  plates  constant  heat  flux  flow 


The  nusselt  number  is  independent  from  the  Peclet  number  for  this  case,  and 

approaches  to  8.235  where  x  goes  to  infinity. 

C.    HEAT  ABSORPTION  UPTO  THE  ENTRANCE 

The  energy  change  from  the  entrance  to  x=0  where  the  heating  starts  has  been 
investigated  and  plotted  for  fully  developed  laminar  flow  inside  the  circular  tube  and  parallel 
plates  for  constant  wall  temperature  and  heat  flux  cases  in  figures  16-19.  The  Peclet  numbers 
0.5,  1,  2,  3,  5,  10,  50,  and  100  were  used  for  the  plots  to  see  how  the  Peclet  number,  Pe=RePr, 
affects  the  heat  absorption  of  the  fluid.  As  it  is  seen  from  the  graphs,  the  heat  absorption 
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logarithmically  increases  as  the  Peclet  number  decreases.  Therefore  axial  conduction  effects  or 
heat  absorption  upto  the  entrance  where  the  heating  starts  has  to  be  considered  for  low  Peclet 
numbers. 


1.    Flow  Inside  the  Circular  Tube,  Constant  Wall  Temperature 


Case 


The  graph  of  the  heat  absorption,  E'  vs.  Peclet  number  for  the  circular  tube, 
constant  wall  temperature  flow  was  plotted  in  figure  16  where  the  dimensionless 
absorbed  energy  E*  is  E*  _      ^absorbed 

{t0-tc)pc/0 

,     ENERGY  ABSORPTION,  CIRCULAR  TUBE.  CONSTANTWALL  TEMP. 
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Figure  16.  The  energy  absoption  upto  the  entrance  for  circular  tube  constant  wall  temperature  flow 

2.    Flow  Inside  the  Circular  Tube,  Constant  Heat  Flux  Case 

The  graph  of  the  heat  absorption,  E*  vs.  Peclet  number  for  the  circular  tube, 
constant  heat  flux  flow  was  plotted  in  figure  1 7  where  the  dimensionless  absorbed  energy 


E*  is  E'  =  — - 


' absorbed 


QoPCprJk 
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ENERGY  ABSORPTION.  CIRCULAR  TUBE,  CONSTANT  HEAT  FLUX 
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Figure  17.  The  energy  absoption  upto  the  entrance  for  circular  tube  constant  heat  flux  flow 


3.    Flow  Between  the  Parallel  Plates,  Constant  Wall  Temperature 


Case 


The  graph  of  the  heat  absorption,  E    vs.  Peclet  number  for  the  parallel  plates,  constant 


wall  temperature  flow  was  plotted  in  figure  18  where  the  dimensionless  absorbed  energy  E    is 


r  = 


' absorbed 


h2bpc(t0-te)' 
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Figure  18.  The  energy  absoption  upto  the  entrance  for  parallel  plates  constant  wall  temperature  flow 
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4.    Flow  Between  the  Parallel  Plates,  Constant  Heat  Flux  Case 

The  graph  of  the  heat  absorption,  E    vs.  Peclet  number  for  the  parallel  plates,  constant 
heat  flux  flow  was  plotted  in  figure   19  where  the  dimensionless  absorbed  energy  E*  is 

E"  *   absorbed 


q0bh  /a 


ENERGY  ABSORPTION.  PARALLEL  PLATES.  CONSTANT  HEAT  FLUX 
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Figure  19.  The  energy  absorption  upto  the  entrance  for  parallel  plates  constant  heat  flux  flow 

D.    AXIAL  CONDUCTION  EFFECTS  INCLUDING  THE  WALL 
CONDUCTION 

The  axial  conduction  effects  for  circular  tube  constant  wall  temperature  case  has 
been  investigated  including  the  wall  conduction  effects.  The  Peclet  numbers  1,  3,  5, 
outer/inner  radius  ratios  of  1.2,  2  and  fluid/wall  thermal  conductivity  ratios  of  1  and  1/10 
were  used  to  see  how  these  properties  affect  the  temperature,  nusselt  number  and  heat 
absorption.  The  theta  bulk  mean  temperature,  nusselt  number  variations  were  plotted  vs. 
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x+  and  the  heat  absorption  were  quantified  as  a  function  of  Peclet  number,  thermal 
conductivity  and  radius  ratios.  The  results  are  seen  on  the  following  pages. 
1.    Bulk  Mean  Temperature  Variations 

The     eci      vs.  x+  graphs  are  seen  in  figures  20. a  and  20. b  for  different  Peclet 

#c/,*=0 

numbers,  thermal  conductivity  and  wall  thickness.  Since  the  bulk  mean  temperature 
values  are  very  small  about  x+=0,  the  normalized  bulk  mean  temperature  on  the  centerline 
is  a  good  way  to  see  the  change  at  the  bulk  temperature  relative  to  the  centerline 

+  n               ,     ,    „                       •     *      V  ~0        ,    +  •             2x  /  d 
temperature  at  x  =0,  where  the  bulk  temperature  is  6  =  -r ^  and  x   is  x    = . 

(t9-h)  RePr 

The  bulk  mean  temperature  vs.  x+  graphs  for  the  Peclet  numbers  1,3,  and  5,  and  the 
thermal  conductivity  ratios  of  1,  1/10  and  the  radius  ratios  of  1.2  and  2  are  seen  in  figure 

(tm-to) 


2 1  where  the  bulk  mean  temperature  is  6m  = 


('.-'J- 


THE  NORMALIZED  BULK  MEAN  TEMPERATURE 


Pe=1  k=1     h=1.2 
Pe=1  k=1/10h=1.2 
Pe=3k=1  h=1.2 
Pe=3k=1/10h=1.2 
Pe=5k=1     h=1.2 
Pe=5k=1/10h=1.2 
Pe=1  without  wall 
Pe=3  without  wall 
Pe=5  without  wall 


Figure  20.a  The  normalized  bulk  mean  temperature  with  and  without  wall 
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THE  NORMALIZED  BULK  MEAN  TEMPERATURE 


Pe=1  k=1     h=2 
Pe=3k=1  h=2 
Pe=5k=1  h=2 
Pe=1  without  wall 
Pe=3  without  wall 
Pe=5  without  wall 


Figure  20.b  The  normalized  energy  absorption  upto  the  entrance  for  parallel  plates  constant  heat  flux  flow 

The  theta  bulk  mean  values  with  the  wall  conduction  case  is  closer  to  the  steady 
state  values  than  the  no  wall  conduction  case  as  it  is  seen  in  figure  2 1 .  The  bulk  mean 
temperature  go  to  zero  as  x  goes  to  plus  infinity.  This  means  that  the  temperature  of  the 
fluid  including  the  wall  conduction  increases  faster  than  the  temperature  of  the  fluid 
neglecting  the  wall  conduction.  The  main  reason  of  the  relative  increase  at  the 
temperature  including  the  wall  conduction  case  is  the  heat  that  is  transferred  by  the  axial 
conduction  of  the  wall  opposite  to  the  flow  direction.  The  heat  moves  to  the  opposite  of 
the  flow  direction  and  is  absorbed  by  the  fluid.  So  the  fluid  temperature  increases  much 
more  than  no  wall  case. 


THE  BULK  MEAN  TEMPERATURE  VS  X 


Figure  21.  The  bulk  mean  temperature  with  and  without  wall  conduction 


66 


2.  Nusselt  Number  vs.  X 

The  Nusselt  numbers  vs.  x+  graphs  are  seen  in  figure  22  for  different  Peclet  numbers, 

thermal  conductivity  and  wall  thickness  where  x+  is  x+  = .  The  Nusselt  number 

RePr 

distribution  is  very  different  from  the  no  wall  conduction  case  as  it  is  seen  in  the  figure. 
As  the  wall  gets  thinner  and  the  thermal  conductivity  increases,  the  values  of  the  Nusselt 
number  gets  closer  to  the  case  of  no  wall  case.  h=l/10  and  k=1.2  is  close  to  the  no  wall 
conduction  case,  so  the  values  of  the  Nusselt  number  for  these  curves  are  close  to  the 
values  of  no  wall  case. 


THE  NUSSELT  NUMBER  VS.  X 
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Figure  22.  The  nusselt  number  distribution  including  the  wall  conduction 
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3.  Heat  Absorption  vs.  X 

The  graphs  of  the  heat  absorption,  E    vs.  Peclet  numbers  for  the  conductivity  ratios  of 
1,1/10  and  radius  ratios  of  1.2  and  2  were  plotted  in  figure  23  where  the  dimensionless  absorbed 


energy  E    is  E*  _      E absorbed      .  As  it  is  seen  from  the  graph,  As  the  wall's  conductivity  or 

thickness  increases,  the  heat  absorption  before  x=0  significantly  increases.  The  wall 
conduction  case  has  to  be  added  to  the  problem  for  the  microheat  exchanger  problems. 
Because  the  relative  thickness  of  the  wall  is  very  large  as  it  is  compared  to  the  inner 
diameter  of  the  tube.  The  values  of  h=1.2  or  h=2  are  realistic  for  the  microheat  machines, 
because  the  diameter  of  the  tubes  are  very  small  for  these  heat  exchanger.  The  heat 
absorption  for  kf/ks=l/10  and  h=1.2  case  is  more  than  ten  times  of  the  case  of  no  wall.  If 
the  wall  thickness  would  be  h=2,  then  the  heat  absorption  would  be  much  greater. 


ENERGY  ABSORPTION  VS.  X 
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k=1  h=2 

-♦_       k=1  h=1.2 

k=1/10h=1.2 
No  wall  case 


10 
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Peclet  Number 


Figure  23.  The  energy  absorption  upto  the  entrance  with  and  without  the  axial  wall  conduction 
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VI.    CONCLUSIONS 


Axial  conduction  effects  in  duct  flows  are  important  and  have  to  be  considered  for 
cases  of  Pe<10.  The  theta  bulk  mean  temperature  and  Nusselt  number  variations  vs. 
downstream  distance  were  plotted  as  seen  in  the  results.  It  is  obviously  seen  the  increase 
in  the  temperature  against  the  flow  direction  from  x=0  as  the  Peclet  number  decreases. 
Heat  absorption  by  the  fluid  in  the  insulated  region  was  also  quantified  by  using  the  non- 
dimensional  parameter  E  .  The  heat  absorption  increases  logarithmically  as  the  Peclet 
number  decreases,  as  seen  in  the  heat  absorption  graphs.  When  the  wall  conduction 
effects  are  included  in  the  problem,  it  was  calculated  that  most  of  the  heat  absorption  by 
the  fluid  is  because  of  the  wall  conductivity  as  the  wall  conductivity  gets  larger.  The  same 
effect  is  also  seen  with  increasing  wall  thickness.  The  heat  absorption  by  the  fluid  gets 
larger  when  the  wall  gets  thicker.  It  is  very  important  to  include  the  wall  conduction  in 
the  problem  as  the  wall  gets  thicker  and  the  thermal  conductivity  of  the  wall  gets  larger. 
With  large  wall  conductivity  and  wall  thickness,  the  fluid  temperature  increases 
significantly  before  the  heating  starts.  The  current  interests  in  axial  conduction  effects  are 
because  of  the  advent  of  micro  heat  exchanger  technology  with  microtubes  where  the 
Reynolds  number  is  very  small. 
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APPENDIX  -  A  SAMPLE  PROGRAM 
The  Program  for  Circular  Tube,  Constant  Wall  Temperature  Case 

This  program  solves  the  finite  difference  equations  for  the  circular  tube,  constant 
wall  temperature  case  by  using  the  Gauss  Seidel  Iteration  Method  with  the  overrelaxation 
parameter  w=1.8.  The  program  solves  the  equations  doing  the  following  steps  : 

1 .  The  matrix  is  formed  by  the  program,  the  first  row  represents  the  region  which 
touches  the  tube.  The  last  row  is  the  centerline.  The  first  column  is  the  condition,  where  x 
goes  to  minus  infinity.  The  last  column  represents  the  boundary  condition  as  x  goes  to 
plus  infinity. 

2.  The  boundary  conditions  are  applied  to  the  problem.  The  centerline  boundary 
condition,  the  boundary  conditions  at  r=r0,  the  boundary  conditions  as  x  goes  to  minus 
and  plus  infinity  are  applied  to  the  last  row,  first  row,  first  column  and  the  last  column 
respectively. 

3.  The  iterations  are  made  calculating  the  values  at  all  the  grid  points.  As  it  is  seen 
in  the  program,  the  odd  rows  are  calculated  first,  the  even  rows  later.  The  convergence  of 
the  problem  is  better  in  this  case. 

4.  The  difference  between  the  error  and  the  boundary  conditions  for  the  plus  and 
minus  infinity  is  checked.  If  the  boundary  conditions  and  the  difference  between  the 
iterations  are  reasonable,  the  program  is  stopped  and  the  data  is  processed.  The  max. 
difference  between  the  iterations  is  less  than  1E-6.  The  max  difference  between  the 
values  of  the  horizontal  grid  points  which  represents  minus  and  plus  infinity  is  less  than 
1E-6. 
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This  problem  solves  the  problem  for  circular  tube 
constant  wall  temperature  case 

Programmed  by  Ibrahim  Girgin 


format  long 

clear 
start =1; 

iteration_number=3  00  000; 

%  Nx  is  the  column  number  of  the  matrix 

%  Ny  is  the  row  number  of  the  matrix 

%  x  zero  is  the  column  number  which  represents  x=0 

Nx=2  501 

Ny=101 

x_zero=1501 

%  For  Peclet  number  is  1 

Pe=l 

%  For  the  relaxation  parameter  1.8 

w=1.8 
save_iteration=100 

hy  is  the  vertical  distance  between  the  grids 
hx  is  the  horizontal  distance  between  the  grids 

hy=l/(Ny-l) ; 
hx=hy ; 

teta  is  the  matrix  which  holds  the  grid  points 
initialization  of  the  teta 

teta=ones (Ny,Nx) ; 

Apply  the  boundary  condition  for  r=Ro,  x>=0 

teta ( : , x_zero :Nx) =zeros (Ny, length (x_zero :Nx) ) ; 


%   e.g. 

J=Ny 

> 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

%  the  initial 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

%  value  of 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

%  theta 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

% 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

% 

j=l 

> 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

i=l  i=x_zero  i=Nx 

The  values  for  this  theta  matrix  is  : 
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%     Nx=14    Ny=6   x_zero=8    hx=0.2    hy=0.2 

%    Pe, w, f ile_name, save  iteration  is  what  is  to  be  given 


%  The  variables  that  are  being  used  inside  the  loop 

varl=PeA2*hxA2; 
var2=2/hyA2, 
var3=l/2/hx; 
var4=l/hy"2 
var5=l/varl-var3 ; 
var6=l/varl+var3 ; 
var7=l/2/hy; 
var8=2/varl+4*var4 ; 
var9=2/varl+var2 ; 
varl0=2*hx; 

%   The  iteration  starts  here.  For  every  value  of  kkk,  one 
%  iteration  is  made 

for  kkk=start : iteration  number 


f or_max_err=teta ; 

for  i=2:Nx-l      %    do  the  iteration  for  all  the  columns  except  the 
first 

%  and  the  last  columns 

for  j=l:2:Ny     %   do  the  iteration  for  all  the  odd  rows 

y= ( j -1) *hy;    %    calculate  the  y,     0<=  y  <=1 

if  j==l       %   do  the  iteration  for  radius=0  (for  the  center) 

teta(Ny,i)= (teta(Ny, i+1) * (var5) +teta (Ny, i-1) * (var6) +  . . . 
4*teta(Ny-l, i) *var4) /var8 ; 

elseif  j==Ny      %    do  the  iteration  for  radius=R 

if  i<x_zero    %    if  i<x_zero  &  for  radius  R,  it  is  insulated 

teta (1, i) =  (  (teta (1, i  +  1) +teta  (1, i- 
1) ) / (varl) +teta (2, i) *var2) /var9; 


end 

else  %  do  the  iteration  for  all  the  odd  rows  except  r=0 

and  r=R 


teta (Ny- j+1, i) =teta (Ny- j+1, i) +w* ( (teta (Ny- j+1, i+1) * (var5  + 
y^2/varl0) . . . 
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+teta (Ny-j+1, i-1) * (var6  -  y^2/varl0) +teta (Ny- j , i) *  (var4  + 
var7/y)  +  .  .  . 

teta (Ny-j+2, i) * (var4  -  var7/y) ) /var9-teta (Ny-j+1 , i) ) ; 


end 

end 
for  j=2:2: (Ny-1)       %     do  the  iteration  for  all  the  even  rows 

y=(j-l)*hy; 

teta (Ny-j+1, i) =teta (Ny-j+1, i) +w* ( (teta (Ny-j+1, i+1) * (var5  + 
y"2/varl0)  .  .  . 

+teta (Ny-j+1, i-1) * (var6  -  yA2/varlO) +teta (Ny- j , i) * (var4  +  var7/y) + 

teta (Ny-j+2, i) * (var4  -  var7/y) ) /var9-teta (Ny-j+1,  i)  )  ; 

end 

end    %     End  of  the  iteration 


%     Save  the  variables  at  every  "save_iteration" 

%     e.g.  if  the  save_iteration  equals  to  20,  then  save 

%         the  variables  at  when  kkk=20, 40, 60, 80 , . . .  so  on. 


if  fix (kkk/save_iteration) ==kkk/save_iteration; 

max_error= (max(max(abs (teta-for_max_err) )  )  ) 
clear  for_max_err 

save  forpel  teta  Nx  Ny  x_zero  w  Pe  hx  hy  max_error  kkk 
save  iteration  start  y 


end 
end 
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